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SUWARY 


A  theory  is  presented  to  calculate  the  flow  field  due  to  main 
and  tail  rotors  of  a  helicopter  in  hovering  flight.  To  solve  this 
complex  flow  problem,  simple  vortex  rings  are  used  to  model  main  and 
tail  rotor  wakes.  Preliminary  numerical  results  are  found  to  be 
encouraging. 
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1.  INTRODUCTION 


In  certain  helicopter  flight  conditions,  the  fuselage, 
and  other  components  such  as  the  tail  rotor,  fin  and  hub,  are 
partially  or  wholly  immersed  in  the  main  rotor  slipstream.  This 
leads  to  interactional  aerodynamic  phenomena  which  have  been  observed 
in  many  helicopters  [1].  A  clear  knowledge  of  these  phenomena  is 
important  not  only  because  they  could  be  a  serious  source  of  vibration 
and  fatigue  [2],  but  also  because  they  may  affect  the  aerodynamic 
efficiency  of  various  components.  The  present  trend  towards  greater 
compactness  necessitates  the  various  aerodynamic  elements  being  close 
to  each  other.  This  means  that  the  aerodynamic  interaction  problem 
will  become  even  more  acute  in  the  future. 

Because  of  the  large  number  of  helicopter  aerodynamic 
components  involved,  a  comprehensive  mathematical  analysis  of  the 
interaction  phenomena  is  extremely  difficult.  Mutual  interaction  of 
the  vortex  elements  generated  by  the  separate  components  will  deform 
the  vortices  which,  because  of  viscosity,  will  eventually  dissipate. 

To  solve  this  complex  interactive  flow  problem,  various  simplifying 
assumptions  are  necessary.  In  particular,  these  include  the  number 
and  geometry  of  the  aerodynamic  elements,  and  the  flow  field  generated 
by  each  element.  In  this  document  only  the  two  important  aerodynamic 
elements,  the  main  rotor  and  tail  rotor,  are  considered. 

The  wake  model  to  represent  a  helicopter  rotor  blade 
varies  from  very  complex  lifting  line  theory  to  simple  blade  tip 
vortex  theory.  As  an  initial  step,  a  simple  tip  vortex  wake  model  is 
used  here  to  represent  the  wake  of  a  rotor  blade.  To  simplify  the 
numerical  process  further,  the  trailing  vortex  pattern  is  assumed  to 
form  a  series  of  vortex  rings.  Similar  wake  models  have  been  used 
previously  by  Cook  [3]  to  calculate  the  helicopter  main  rotor 
performance,  and  recently  Rayner  [4]  used  a  ring  vortex  wake  model  to 
study  hovering  flight  of  animals.  Other  forms  of  representing  the 
wake  such  as  the  skewed  spiral,  necessitate  numerical  integration  and 
are  therefore  time  consuming. 

In  developing  the  theory,  consideration  has  been  limited 
to  aircraft  having  two  rotors,  although  no  restriction  has  been 
placed  on  the  number  of  rotor  blades.  A  simple  and  efficient 
procedure,  suitable  for  flight  simulation  studies,  for  the  calculation 
of  the  flow  field  due  to  main  and  tail  rotors  is  described.  For  the 
test  case,  the  induced  flow  in  the  plane  of  the  main  rotor  in  steady 
hovering  flight  has  been  calculated  and  compared  with  results  from  a 
more  complex  model.  The  theory  can  be  readily  extended  to  cope  with 
forward  flight.  Throughout,  it  is  assumed  that  the  fluid  is  inviscid 
and  incompressible. 

2.  DEVELOPMENT  OF  THE  WAKE  MODEL 

To  estimate  the  combined  main  and  tail  rotor  velocity 
field,  the  wake  of  each  rotor  is  modelled  as  a  stack  of  circular 
coaxial  vortex  rings  beneath  the  rotor.  The  rings  are  equally  spaced, 


the  first  one  being  half  the  spacing  distance  below  the  rotor.  The 
momentum  theory  is  used  to  calculate  the  spacing  between  the  rings. 


Initially,  all  the  rings  were  the  same  radius  as  the 
rotor  so  as  to  simplify  the  situation  but  a  comparison  with  results 
from  other  models  (Reference  5)  showed  a  large  variation  in  results. 
The  radii  of  the  rings  are  now  given  by  - 


R  *  (0.78  +  0.22  e"**)  R* 


(see  Reference  6)  where  K  is  a  wake  contraction  parameter  (*4(CT)*), 
-  2*0-1)  for  the  ith  vortex  ring  and  R'  is  the  radius  of  the 
rotor.  This  gives  much  improved  results  for  induced  velocity  in  the 
plane  of  the  rotor,  except  for  r/R1  >  0.975  -  the  very  edge  of  the 
rotor  disc. 


The  vortex  rings  are  assumed  to  be  sufficiently  far 
apart  so  that  they  have  no  effect  on  one  another.  Thus,  the  total 
velocity  induced  by  the  entire  "stack"  of  rings  is  just  the  vector 
sum  of  the  velocities  Induced  by  the  individual  rings. 

The  tail  rotor  is  modelled  in  exactly  the  same  way  as 
the  main  rotor,  although  a  different  co-ordinate  system  is  used,  and 
the  results  can  then  be  combined  with  those  from  the  main  rotor  to 
give  the  total  induced  velocity.  For  the  data  used,  however,  the 
contribution  to  the  induced  velocity  from  the  tail  rotor  is  small 
except  near  the  rear  edge  of  the  main  rotor  disc. 

3.  DESCRIPTION  OF  THE  COMPUTER  PROGRAM 


The  program  itself  is  a  straight-forward  implementation 
of  the  formulae  developed  in  the  Appendix,  The  required  inputs  are: 

the  radius  of  the  rotor 

the  separation  between  consecutive  vortex  rings 
the  number  of  vortex  rings  to  consider 
the  circulation  of  each  ring. 

If  the  effects  of  a  tail  rotor  are  to  be  included  then  this  informa¬ 
tion  is  also  needed  for  the  tail  rotor,  as  well  as  the  position  of 
the  centre  of  the  tail  rotor  relative  to  the  main  rotor. 

The  Induced  velocity  at  the  rotor  is  then  found  using 
the  information  provided,  by  application  of  the  formulae  in  the 
Appendix.  If  the  tail  rotor  is  Included,  the  Induced  velocity  due 
to  each  stack  of  rings  is  found  and  the  results  combined  to  give  the 
total  effect.  In  this  case,  the  effect  of  each  rotor  is  output,  as 
well  as  the  total  effect. 


The  units  of  the  output  are  determined  by  those  of  the 
input  data.  If  distances  (i.e.  radii  and  separations)  are  given  in 
metres  and  circulations  in  metres2/ second,  then  the  results  are  in 
m/s,  but  if  other  output  units  are  desired  then  the  input  data  must 
be  modified  appropriately. 

The  induced  velocities  are  calculated  using  two  right- 
handed  Cartesian  coordinate  systems.  The  main  one,  to  which  all 
input  and  output  is  referenced,  is  centred  at  the  centre  of  the  main 
rotor,  has  the  x-axis  directed  horizontally  backwards,  the  y-axis 
horizontal  and  the  z-axis  vertically  upwards.  The  second  system, 
which  is  used  for  convenience  in  calculating  induced  velocities  due 
to  the  tail  rotor,  is  parallel  to  but  rotated  from  the  main  system 
as  shown  in  Figure  2. 

If  the  centre  of  the  second  co-ordinate  system  is 
(a,b,c)  relative  to  the  main  system,  then  the  point  P(x,y,z)  in  the 
main  system  is  given  by  P'(x-a,z-c,  b-y)  relative  to  the  second 
system. 

The  program  was  written  in  modular  form,  so  that  the 
main  program  controls  the  flow  of  execution  but  does  very  little  of 
the  actual  work;  most  of  the  work  is  done  in  subroutines  to  read 
in  the  input  values  and  calculate  the  induced  velocities  by  applica¬ 
tion  of  the  formulae.  Note  that  the  calculations  require  the  use  of 
a  numerical  integration  routine  from  the  library  DSKB.NAGF  [7,7], 
which  must  be  included  at  link  time  for  the  program  to  run  correctly. 

4.  RESULTS  AND  CONCLUDING  REMARKS 

The  wake  in  the  hovering  flight  of  a  rotor  is  modelled 
by  a  stack  of  circular  blade  tip  vortex  rings.  Experimental 
evidence  suggests  that  blade  tip  vortex  is  the  most  dominent  element 
in  the  rotor  wake.  Experimental  results  also  show  that  this  blade 
tip  vortex  though  highly  distorted  farther  away  from  the  rotor, 
follows  a  path  quite  close  to  a  circle  in  the  immediate  vicinity  of 
the  rotor.  If  calculation  of  the  velocity  field  close  to  the 
rotor  is  of  main  concern  then  representation  of  the  rotor  vortex 
wake  by  a  stack  of  vortex  rings  is  physically  reasonable.  This  also 
simplifies  the  mathematical  analysis,  the  computer  time  is  reduced 
which  make  the  model  very  suitable  for  flight  simulation  studies. 

Induced  velocities  in  the  rotor  plane  calculated  using 
the  ring  vortex  wake  model  are  presented  in  Figure  3  for  a  typical 
four  bladed  helicopter  rotor  in  hover.  A  comparison  of  results  from 
this  model  with  those  in  Reference  5  shows  reasonable  agreement 
between  the  theories  except  at  the  edge  of  the  blade,  where  some  of 
the  assumptions  made  deriving  the  current  formulae  may  not  hold. 
Methods  to  improve  the  present  mathematical  model  are  being  investi¬ 
gated.  The  theory  can  be  readily  extended  to  cover  forward  flight 
of  any  two  rotor  helicopter. 


REFERENCES 


P.F.  Sheridan  and  R.P.  Smith,  Interactional  Aerodynamics  - 
A  New  Challenge  to  Helicopter  Technology,  Journal  of 
American  Helicopter  Society,  January  1980. 

J.W.  Leverton,  J.S.  Pollard  and  C.R.  Wills,  Main  Rotor  Wake/ 
Tail  Rotor  Interaction,  VERTICA,  Vol.  1,  1977,  pp  213-221. 

C.V.  Cook,  Induced  Flow  Through  a  Helicopter  Rotor  in  Forward 
Flight,  Research  Paper  374,  Westland  Helicopters  Ltd., 
Somerset,  1970. 

J. M.V.  Rayner,  A  Vortex  Theory  of  Animal  Flight.  Part  1. 

The  Vortex  Wake  of  a  Hovering  Animal,  J.  Fluid  Mech.,  Vol.  91 
1979,  pp  697-730. 

K. R.  Reddy,  Prediction  of  Helicopter  Rotor  Downwash  in  Hover 
and  Vertical  Flight,  ARL  Aerodynamics  Report  150,  January 
1979. 

J.D.  Kocurek  and  J.L.  Tangier,  A  Prescribed  Wake  Lifting 
Surface  Hover  Performance  Analysis,  American  Helicopter 
Society  32  V/STOL  Forum,  May  1976. 


APPENDIX 


MATHEMATICAL  DERIVATION  OF  THE  VELOCITY  FIELD  INDUCED  BY 
A  VORTEX  RING  ELEMENT  (from  reference  4) 


Consider  the  effect  of  a  single  vortex  ring,  centred  at  the 
origin  and  lying  in  the  x-y  plane,  on  a  point  in  the  x-z  plane 
(fig.  la).  No  generality  is  lost  by  this  positioning  of  the  point 
as  any  point  can  be  so  located  by  a  rotation  of  the  co-ordinate 
system  about  the  z-axis. 


As  long  as  the  point  is  not  near  the  core  of  a  vortex  ring 
the  vortex  core  can  be  approximated  by  a  line  vortex,  and  the  induced 
velocity  field  is  given  by  the  Biot-Savart  law. 


S(i) 


dL  (x-X(s))At(s) 
4*  i |x-X(s)| Is 


(1) 


where  s  is  a  parameter  of  the  circumference  of  the  vortex  ring. 


In  the  case  of  figure  la  with  points  P  and  Q  (assumed  to  be 
in  the  x-z  plane). 


X  =  R(cose,  sine,  0) 

x  -  (r,o,h)  (2) 

s  *  Re 

t  «  M  *  (-sine,  cose,  0) 
ds 


so  from  (1)  it  follows  that 

2w  (hcose,  hsine,  R-rcose) 

[rz-2rRcose  +  R2  +  h2]3/3 
o 


(3) 


a  a  a 

Denote  the  unit  vectors  in  the  x,y  and  z  directions  by  x,y,z 
respectively. 


By  symmetry  about  the  z-axis 


Define 


u.y  *  0 


u.x  *  ^  I2(r,R,h) 


u.z  =  (RIi (r,R,h)  -  rI2(r,R,h)) 


where  Ij  and  I2  can  be  evaluated  in  terms  of  the  complete  elliptic 
integrals 

K(e)  =  [  ^2(l-e2cos2e)"*  de 


E(e)  *  [  2(l-e2cos2e)*  de 


where  the  eccentricity  e  is  given  by 
e2  =  4rR[(r+R)2  +  h2]'1 


For  K(e)  to  be  real,  e2  <  1  which  is  always  true  as  follows 


(r-R)2  +  h2  >  0 


V  r,R,h 


r2  +  2rR  +  R2  +  h2  >  4rR 
(r+R)2+  ha  a  4rR 


[(r+R)2  +  h2] 


e2  i  1 


Comparing  (3)  and  (5)  gives 

-  *  kRh  f2ir  cose  de 


_  [(r+R)2+h2-2rR-2rRcose]3 


*  Ig  *  cose[(r+R)2+h2-2rR(l+cose)]  de 


•3/2  p 

I2  a  C(r+R)2+h2]  *  cose  |i.e^lll£2se)]~3/2  d0  (10) 

.  o 


But  cos2e  * 

2 


17  *  2[(r+R)2+h2]‘3/2  P  cos2e[l-e2cos2e]'3/2  de  (11) 


Also  note 

[  cos2e[l-e2cos2e]  ^2  de  =  2 (  cos2e[l-e2cos2e]  ^  de 


I,  *  4[(r+R)2+h2]’3^2  f  ^2cos2e[l-e2cos2e]“3/r2  de 


4H(e)[(r+R)2+h2] 


2+h21"3/2 


from  tables  of  Integrals,  where  H(e)  is  defined  by 


H(e)  «  1-  f  2^4  E(e)  -  2K(e)l 
e2  ll-e2  J 


Substituting  (9)  back  into  (5)  and  comparing  the  result  to  (3) 
and  proceeding  similarly  to  the  above  gives 


Ij  »  4[(r+R)l+h2]"3^2  ^2(l-e2cos2e)"3^2  de 

•  o 

»  46(e)  [(r+R)2+h2r3/2 


(18) 


where  G(e)  = 


(1-e2) 


From  (14)  and  (17)  it  can  be  seen  that 


u  =  u  x  +  u  z 
r  z 


kRh 


r  4n  [(r+R)2+h2]3/2 


H(e) 


kR 


z  4ir  [(r+R)2+h2]3/2 


(RG(e)-rH(e)) 


(19) 

(20) 

(21) 


with  e  given  as  in  (7). 


If  the  point  Q  does  not  lie  in  the  x-z  plane  then  the  quantity 
ur  is  a  radial  velocity  directed  outwards  from  the  z-axis,  which 
must  be  split  into  true  x  and  y  components. 

Then 


u  x  *  ur  cos  4 
u  y  =  ur  sin  $ 


(22) 


where  4  is  the  angle  shown  in  figure  lb. 
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